Roughening Transition in a Moving Contact Line 
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The dynamics of the deformations of a moving contact line on a disordered substrate is formulated, 
taking into account both local and hydrodynamic dissipation mechanisms. It is shown that both the 
coating transition in contact lines receding at relatively high velocities, and the pinning transition 
for slowly moving contact lines, can be understood in a unified framework as roughening transitions 
in the contact line. We propose a phase diagram for the system in which the phase boundaries 
corresponding to the coating transition and the pinning transition meet at a junction point, and 
suggest that for sufficiently strong disorder a receding contact line will leave a Landau-Levich film 
immediately after depinning. This effect may be relevant to a recent experimental observation in a 
liquid Helium contact line on a Cesium substrate [C. Guthmann, R. Gombrowicz, V. Repain, and 
E. Rolley, Phys. Rev. Lett. 80, 2865 (1998)]. 



I. INTRODUCTION AN SUMMARY 

When a drop of liquid spreads on a solid surface, the 
contact line, which is the common borderline between 
the solid, the liquid, and the corresponding equilibrium 
vapor, undergoes a rather complex dynamical behavior. 
This dynamics is determined by a subtle competition 
between the mutual interfacial energetics of the three 
phases, dissipation and hydrodynamic flows in the liq- 
uid, and the geometrical or chemical irregularities of the 
solid surface jjj. 

For a partially wetting fluid on sufficiently smooth 
substrates, a contact line at equilibrium has a well de- 
fined contact-angle 9 e that is determined by the solid- 
vapor 7sv and the solid-liquid 7sl interfacial energies, 
and the liquid surface tension 7 through Young's rela- 
tion: 7sv — 7sl = 7COs6' e . For a moving contact line, 
however, the value of the so-called dynamic contact-angle 
9 changes as a function of velocity: 8 > 9 e for an advanc- 
ing contact line and 9 < 9 e for a receding one. This is be- 
cause the unbalanced interfacial force 7sv — 7sl — 7 cos 9 
now has to be balanced with a frictional force in a steady 
state situation. The dissipation in the moving contact 
line, which is responsible for the friction, can be either 
of local origin, for example due to microscopic jumps of 
single molecules (from the liquid into the vapor) in the 
immediate vicinity of the contact line [^],[| , or due to vis- 
cous hydrodynamic losses inside the moving liquid wedge 

For a contact line that is receding at a velocity v, it 
has been shown by de Gennes |J that a steady state is 
achieved in which the liquid will partially wet the plate 
with a nonvanishing dynamic contact-angle 9 only for ve- 
locities less than a certain critical value. The dynamic 
contact- angle decreases with increasing v, until at the 
critical velocity the system undergoes a dynamical phase 



transition in which a macroscopic Landau-Levich liquid 
film [^|J^] , formally corresponding to a vanishing 9, will re- 
main on the plate. One can think of the dynamic contact- 
angle as the order parameter characterizing this coating 
transition, in analogy with equilibrium phase transitions, 
while velocity is playing the role of the tuning parame- 
ter. Elaborating further on this analogy then seems to 
suggest that the nature of the coating transition depends 
crucially on the dominant dissipation mechanism: In the 
local picture 9 vanishes continuously as v approaches the 
critical velocity, which makes it look like a second order 
phase transition, while on the contrary, in the hydrody- 
namic picture a jump is predicted in 9 from 9 e /y/3 to 
zero at the transition, which is the signature of a first 
order phase transition Q . 

Another notable feature of contact lines, which is re- 
sponsible for their novel dynamics, is their anomalous 
long-ranged elasticity jnj. For length scales below the 
capillary length (which is of the order of 3 mm for wa- 
ter at room temperature), a contact line deformation of 
wavevector k will distort the surface of the liquid over 
a distance Assuming that the surface deforms in- 

stantaneously in response to the contact line, the elastic 
energy cost for the deformation can be calculated from 
the surface tension energy stored in the distorted area, 
and is thus proportional to |fc|. The anomalous elasticity 
leads to interesting equilibrium dynamics, corresponding 
to when the contact line is perturbed from its static po- 
sition, as studied by de Gennes Balancing the rate 
of interfacial energy change and the dissipation, which 
he assumed for small contact-angles is dominated by the 
hydrodynamic dissipation in the liquid nearby the con- 
tact line, he finds that each deformation mode relaxes 
to equilibrium with a characteristic (inverse) decay time 
r _1 (fc) = c\k\, where c is a characteristic relaxation ve- 
locity pj . The relaxation is thus characterized by a dy- 
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namic exponent z, defined via r~ 1 (k) ~ |fc| z , which is 
equal to 1. The linear dispersion relation implies that a 
deformation in the contact line will decay and propagate 
at a constant velocity, as opposed to systems with normal 
line tension elasticity, where the decay and the propaga- 
tion are governed by diffusion. This behavior has been 
observed, and the linear dispersion relation has been pre- 
cisely tested, in a recent experiment by Ondarcuhu and 
Veyssie (l2). 

In reality, the presence of defects and heterogeneities 
in the substrate, which could be due to (surface) rough- 
ness or chemical contamination, further complicates the 
dynamics of a contact line ||Tl||l4j| . In the presence of 
such heterogeneities, a contact line at equilibrium be- 
comes rough, because it tries to locally deform so as to 
find the path with optimal pinning energy Jl]] . This is in 
contrast to the case of a perfect solid surface, where the 
contact line is flat. The roughness can be characterized 
as a scaling law that relates the statistical width W of the 
contact line to its length L, via W ~ L^. The so-called 
roughness exponent £ is equal to 1/3 for a contact line 
at equilibrium on a surface with short-range correlated 
disorder |15|-|l8|]. Since the contact line is pinned by the 
defects, a nonzero (critical) force is necessary to set it 
to motion, through a depinning transition ]T9|-p2"| . For a 
contact line at the depinning threshold, a roughness ex- 
ponent about 0.4 has been predicted theoretically from 
two-loop field theoretical renormalization group calcula- 
tions pi! , and numerical simulations |p4| , which seems 
to disagree with the experimental finding of 0.5 for both 
liquid Helium on a Cesium substrate fl25j and water on 
glass experiments pf|. It is also important to note that 
there may be numerous metastable states for the contact 
line due to the random disorder, leading to hysteresis in 
the contact-angle 




FIG. 1. A contact line moving on a disordered substrate 
undergoes shape fluctuations. A portion of the contact line 
can be instantaneously pinned by a defect, thereby nucleating 
domains of typical sizes given by a correlation length £. These 
domains become rough, because they have to conform to the 
minimum energy configuration on the substrate. At the onset 
of a roughening transition this correlation length diverges. 



Here we study the nonequilibrium dynamics of the de- 
formations of a moving contact line on a disordered sub- 
strate p7| , p8fl . The dynamics is governed by a balance 
between three different forces: (i) the interfacial force, 
(ii) the frictional force which can stem from either local 
or hydrodynamic dissipation processes, and (iii) a ran- 
dom force caused by the disorder. We find that the re- 
laxation spectrum of a moving contact line is the same 
as the equilibrium case, but the characteristic relaxation 
velocity depends on v: It decreases with v until at the 
critical velocity corresponding to the coating transition 
it vanishes identically. The progressively slow relaxation 
of a distorted contact line near the coating transition is 
in agreement with a nucleation picture of the phase tran- 
sition (See Fig. |l|). 

We find that coating transition can be actually under- 
stood in terms of a roughening transition of the contact 
line on the disordered substrate. Since linear relaxation 
becomes infinitely slow in the vicinity of the coating tran- 
sition, the dominant relaxation is thus governed by non- 
linear terms in the dynamical equation, and the dynam- 
ical phase transition can thus be properly accounted for 
only by using systematic renormalization group (RG) cal- 
culations. We find that disorder favors the coating tran- 
sition, in the sense that the onset of leaving a Landau- 
Levich film for a random substrate with strength g takes 
place at a dynamic contact-angle 
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for hydrodynamic dissipation, to the leading order. (cv c / 
and a c h are numerical constants to be defined below.) 
The value of the roughness exponent at the transition is 
found to determine the order of the transition. Although 
we find that this exponent acquires non-universal values, 
it appears that the predicted nature of the phase tran- 
sition from the RG calculation is in agreement with the 
mean-field results, i.e. second order for local case and 
first order for hydrodynamic case, for sufficiently weak 
disorder. 

For a sufficiently slowly moving contact line, we find 
that a pinning transition takes place, which is the re- 
verse of the depinning transition, and that it can also be 
understood as a roughening transition. We obtain the 
phase boundary for this transition corresponding to the 
advancing and receding contact angles at the onset of 
pinning as 



el 



i = i 



Or 



9 

101 



(3) 



for the local case, and 
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and some discussions in Sec. IX. Some details related 



for the hydrodynamic case, to the leading order, which 
is in agreement with a previous prediction by Robbins 
and Joanny []l6|| . (a p i and a p h are numerical constants 
to be defined below.) We further find that the roughness 
exponent at the pinning threshold is also non-universal. 

We combine our results for the coating transition and 
the pinning transition, and propose a phase diagram for 
contact lines with local dissipation as depicted in Fig. |^, 
and a corresponding one for contact lines with hydrody- 
namic dissipation as depicted in Fig. |j| In particular, we 
find that the phase boundaries corresponding to the coat- 
ing transition and the pinning transition meet at a junc- 
tion point T, and suggest that for sufficiently strong dis- 
order a receding contact line will leave a Landau-Levich 
film immediately after depinning. This corresponds to 
the dashed lines in Figs, p and ^[ Note that the asymp- 
totic form for the coating transition lines in Figs. and 
H are given in Eqs. (Q) and (||), and the asymptotic form 
for the pinning transition lines are given in Eqs. (§) and 
(||), respectively. 
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FIG. 2. The suggested phase diagram of a contact line with 
local dissipation on a disordered substrate. The asymptotic 
forms for the coating and pinning transition lines are given in 
Eqs. (|l]) and (^), respectively. 

The rest of this article is organized as follows. In Sec. 
H the main ingredients in the dynamics of the contact 
line are discussed, and they are put together in Sec. Ill 
where a stochastic dynamical equation is proposed. In 
Sec. [i"v| , the stochastic dynamical equation is character- 
ized by its self-affine behavior in terms of various expo- 
nents. The mean-field theory of a moving contact line is 
discussed in Sec . M , accompanied by a linear relaxation 
theory in Sec. \J% The effect of the nonlinearities is 
incorporated u sing a dynamical RG scheme which is dis- 
cussed in Sec. VII, followed by the results in Sec. VIII 



to the differential geometry of the moving liquid drop is 
relegated to Appendix [X| and finally some useful asymp- 
totic forms of the results presented in Sec. |l| are given 
in Appendix 0. 
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FIG. 3. The suggested phase diagram of a contact line 
with hydrodynamic dissipation on a disordered substrate. The 
asymptotic forms for the coating and pinning transition lines 
are given in Eqs. (^) and ([|), respectively. Note that the 
coating transition line starts at l/\/3 ~ 0.577 for zero disor- 
der. 



II. DYNAMICS OF A DEFORMING CONTACT 
LINE 



Let us assume that the contact line is oriented along 
the x-axis, and is moving in the y-direction with the po- 
sition described by y(x,t) = vt + h{x,t), as depicted in 
Fig. @. 



A. Interfacial Forces 



If a line element of length dl = dx^J\ + (d x h) 2 is dis- 
placed by 5y(x,t), the interfacial energy will be locally 
modified by two contributions: (i) the difference between 
the solid-vapor 7sv arid the solid-liquid 7sl interfacial 
energies times the swept area in which liquid is replaced 
by vapor, namely, (7sv - !Sh)dlSy/ ^1 + {d x h) 2 , and (ii) 
the work done by the surface tension force, whose di- 
rection is along the unit vector T that is parallel to the 
liquid-vapor interface at the contact and perpendicular 
to the contact line, as 7T -ydlSy. Note that we are inter- 
ested in length scales below the capillary length, where 
gravity does not play a role. The overall change in the 
interfacial energy of the system can thus be written as 
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SE = / dxy/l + (d x h) 2 



7 cos 8 e 



7 T-y 



8y(x,t), 
(5) 



in which we have made use of the Young equation. Note 
that both "forces" should be projected onto the y-axis 
when calculating the work done for a displacement in 
this direction. 



small neighborhood the amount of dissipation is com- 
pletely determined by the local value of the contact line 
velocity, while all the molecular details of the dissipation 
is encoded in an effective friction coefficient The 
overall dissipation can then be written as 



Pi 



1 



dxy/l + [d x hf 



v + d t h(x,t) 



(7) 




FIG. 4. The schematics of the system. 



B. Dissipation 

To calculate the dissipation relevant to the dynamics 
of the contact line, we should consider only the normal 
component of the velocity If the contact-angle is 
very small, the dominant contribution to the dissipation 
comes from the viscous losses in the hydrodynamic flows 
of the liquid wedge [Q. For a slightly deformed contact 
line, we can assume that the dissipation is well approx- 
imated by the sum of contributions from wedge-shaped 
slices with local contact- angles 9(x, t), as shown in Fig. ||. 
This is a reasonable approximation because most of the 
dissipation is taking place in the singular flows near the 
tip of the wedge HJllj]. Using the result for the dissipa- 
tion in a perfect wedge which is based on the lubrication 
approximation p],p9| , we can calculate the dissipation in 
the hydrodynamic regime as 



dx^l + (d x hf 



3rj£ 
6{x,t) 



v + d t h(x, t) 



(6) 



where rj is the viscosity of the liquid and £ is a logarithmic 
factor of order unity ||ll| . One can show that the error in 
the above estimate, which comes from overlap between 
the neighboring slices, only leads to curvature terms that 
are sub-dominant in the long wavelength limit. 

Another physical process that is involved in causing 
dissipation is molecular jumps near the tip of the con- 
tact line, and is local in nature Therefore, in any 



C. Force Balance 

We can find the governing dynamical equation by 
balancing the total friction force obtained as 8 (Pi + 
Ph) I '8d t h(x,t) with the interfacial force —8E/8h(x,t) at 
each point along the contact line. In the limit of small 
contact- angles, we find |30[| 
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0(x,t) 



v + d t h(x,t) 7 
Vl + (cU) 2 ~ 2 



[e 2 e -e(x,tf 



(8) 



The above equation might simply be recovered by locally 
applying the result of Ref. J31J for straight contact lines, 
with the additional geometrical factor (that is needed 
when the direction of motion is not perpendicular to the 
contact line ||) taken into account. 

We can introduce a characteristic velocity for the hy- 
drodynamic friction as Co/, = jOg/^rji), and a corre- 
sponding velocity for the local friction as cqi — /J.j0 2 . It 
is then useful to write the dynamical force balance equa- 
tion in terms of these characteristic velocities. It reads 
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f 



cqi c oh 6(x,t) 



v + d t h(x, t) 
VI + (d x h) 2 



6( Xl tf 



(9) 



Note that the relative strength of the two dissipation 
mechanisms is characterized by the ratio coh/coi, and 
we can readily obtain the asymptotic form of the equa- 
tion when local dissipation is dominant by taking the 
limit Coh/coi —> °Oi an d the corresponding form when 
hydrodynamic dissipation dominates by taking the limit 
coh/coi -> 0. 



D. Solving for the Surface Profile 

To complete the calculation, we need to solve for the 
profile of the surface, and in particular, the angle 9(x,t) 
as a function of h(x, t). 

We may assume that the pressure at the surface equi- 
librates rapidly enough, so that for the effective study of 
the long time dynamics of the contact line it will be suf- 
ficient to set the instantaneous Laplace pressure to zero 
For sufficiently small contact-angles, the surface 
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profile z(x,y,t) near the contact line can then be found 
as a solution of the Laplace equation 

V 2 z(x,y,t) = 0. (10) 

One can write a general solution for Eq. (|l(]) , of the form 

/dk 
—/3(k,t)exp \ikx~ \k\(y-vt)] 

(11) 

that yields the moving liquid wedge profile for a flat con- 
tact line. Imposing the boundary condition z(x, vt + 
h(x,t),t) = at the position of the contact line then 
yields 



Combining the above equations then yields 



/3(k,t) = -e 



h(k,t) + J ^-\q\h(q,t)h(k - q,t) + 0(h 3 



(12) 



Note that we have performed the calculation up to the 
second order in h to find the leading order nonlinearity 
in the dynamical equation, and that since the nonlinear 
terms neglected in the expression for Laplace pressure 
are of the third order in z, this calculation is consistent. 

From the slope of the liquid surface at the position of 
the contact line (see Fig. [|and Appendix [A]) one can then 
obtain an expression for the contact-angle as a function 
of the contact line deformation. We find 



{x, t) = 6 { 1 + b 



dk 
2^ 



\k\h(k,t)e 



■ikx 



2 / 2tt 2tt 1 



+ b 2 \k\\k'\ 
+b 3 \k + k'\(\k\ + \k'\-\k + k'\)} 
x h{k,t)h(k',t)e^ k+k '^ x \ , 



(13) 



with bg = 1, b\ = 1, 62 = 0, and 63 = 1. This can then 
be used in Eq. (^) to yield the dynamical equation. 

One may, however, question the validity of the instan- 
taneous pressure relaxation assumption. To improve on 
this approximation, one should attempt to solve for the 
dynamics of the liquid surface together with the contact 
line dynamics, and examine the corresponding time scales 
for the surface and contact line relaxations. 

This dynamics can be formulated, within the frame- 
work of the lubrication approximation, using a continuity 
equation of the form 



(14) 



d t z(x,y,t) = V • ( — Vp 



where the fluid film is locally described as a Poiseuille 
flow under the influence of the gradient of the Laplace 
pressure 



d t z(x,y,t) 



0, 



(16) 



which is the dynamical equation for the surface defor- 
mation in the lubrication approximation ||. To proceed 
systematically, one should attempt to solve Eq. (|l6|) for 
a moving contact line with equilibrium contact-angle 6 e 
subject to the boundary condition z(x, vt+h(x, t), t) = 0, 
perturbatively in powers of the contact line deformation 
h up to second order. Unfortunately, this seems to be 
a formidable task, because of the complex structure of 
this nonlinear partial differential equation. In fact even 
at the zeroth order, i.e. for the case of a flat contact 
line, this problem is still the subject of much theoretical 
investigation @,|,|,|9]|32 1 . 

We can instead try to estimate the surface relaxation 
time from Eq. ( |l6| ) using dimensional arguments. If we 
consider a deformation of the characteristic size q^ 1 (in 
both x and y directions), and put in a wedge- like profile 
in the nonlinear term of Eq. ( [i"6| ) of the form z ~ Oq^ 1 , 
we find that the corresponding (inverse) relaxation time 
scales as r^ 1 (q) ~ (-f9 3 /ri)q. Interestingly, this is the 
same as the time scale that we find for the contact line 
relaxation [see Eq. fl48| ) below], and it shows that the 
instantaneous surface relaxation assumption is not feasi- 
ble. However, since surface relaxation introduces no new 
time or length scales in the system, it seems plausible 
to assume that the contact-angle profile as a function of 
the contact line deformation, as obtained from a full sys- 
tematic solution of Eq. (|l6|), will still maintain the form 
given in Eq. ( [j"3| ) in the long time and long length-scale 
limit, perhaps with different values for the numerical co- 
efficients b n . Since we will be interested only in this limit 
in the context of the RG calculations, it seems reasonable 
to use the general form proposed in Eq. (p"3|). 



E. Disorder 

In most practical cases, the dynamics of a contact line 
is affected by the defects and heterogeneities in the sub- 
strate, in addition to dissipation and elasticity that we 
have considered so far. If the interfacial energies 7sv and 
7sl are space dependent with the corresponding averages 
being 7sv and 7sl, a displacement Sy(x, t) of the contact 
line is going to lead to a change in energy as 



SEd = / dxg(x,vt + h(x,t))Sy(x,t), 



P(x,y,t) = --/V 2 z(x,y,t). 



(15) 



(17) 
where 

g(x,y) = jsv{x,y) -j S h(x,y) - (7sv -7sl)- (18) 

Incorporating this contribution in the force balance leads 
to an extra force term g(x,vt) on the right hand side of 
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Eq. (g) , which would act as a noise term in the dynamical 
equation for contact line deformation of the form 



T)(x,t) 



(19) 



to the leading order. 

Assuming that the surface disorder has short range cor- 
relations (so that the correlation length is a microscopic 
length a) with a strength g, with a Gaussian distribution 
described by 

(g(x,y)) = o, 

(g(x, y)g(x / , y')) = g 2 a 2 S(x - x')5{y - y'), (20) 
we can deduce the distribution of the noise as 



(r)(x, t)r)(x',t')) = 2DS(x - x')5(t - t'), 



where the strength of the noise is given as 

inr I ntf \ 
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and the nonlinear coupling constants 
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(27) 



(28) 



(29) 



III. DYNAMICAL EQUATION OF MOTION 

We can now put together all the different ingredients 
of the dynamics of the contact line that we discussed in 
Sec. ||, an d obtain the governing dynamical equation. 
Inserting Eq. ( |l3| ) in the force balance relation Eq. (|^), 
and adding the noise term of Eq. (119) , we find 



d t h(k, t) = -c\k\h(k, t) + r){k, t) 

A J S A(9, k ~ q)h{q ' t)h{k ~ q > <} ' (23) 

up to the second order in deformations, with a corre- 
sponding nonlinear coupling as 



X(q, k-q) = -Xiq(k - q) + \ 2 \q\\k - q\ 
+X 3 \k\(\q\ + \k-q\-\k\). 



(24) 



The zeroth order term in the force balance equation pro- 
vides the contact-angle-velocity relation as 



in terms of the dynamic contact-angle. The spectrum of 
the noise term in Fourier space is characterized as 

(n(k,t)) = o, 

(r]{k, t)r](k', t')) = 2D(2Tr)S(k + k')S(t - t'), (30) 

where the strength of the noise D is given as in Eq. , 
or equivalently as 
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(l + £0h±\ 

y 1 + c , 6 e ) 



(31) 



The above dynamical equation and its corresponding 
physical implications will be discussed in detail in the fol- 
lowing sections. Since we will be attempting to compare 
the two alternative dissipation scenarios during these dis- 
cussions, we have summarized in Appendix ^ the limit- 
ing form of the above equations corresponding to each of 
these mechanisms separately. 
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(25) IV. CHARACTERIZING THE STOCHASTIC 

DYNAMICS 



and we can also find the other coupling constants, namely 
the relaxation speed 
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Due to the presence of the heterogeneities in the sub- 
strate, the contact line undergoes dynamical fluctuations 
during its (average) drift motion. These fluctuations, 
which are governed by the stochastic dynamical equation 
given in Sec. |l| [Eq. (|23|)], can best be characterized by 
the width of the contact line, which is defined as 
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W 2 (L,t) 



L" 



d d x(/i(x, t) 2 



d d k 



(2tt) 



i<IMM)l' 



(32) 



where the averaging is with respect to the noise term in 
Eq. ( p3] ) . Note that we have generalized the contact line 
to a d-dimensional object of size L so that the dependence 
on the dimensionality becomes manifest. 

Since the stochastic dynamics described by Eq. ( p3| ) 
corresponds to scale-free fluctuations, the resulting two- 
point correlation function should have a scaling form as 



(IMM)| S 



i 



Q{kH). 



(33) 



The function Q (u) has the property that it saturates to a 
finite value for large it, to ensure that a stationary regime 
can be achieved in the long time limit. Using the above 
scaling form, the width of the contact line will be given 
as 



W 2 (L,t) 



k^^kH), 



tt/L 



(34) 



that yields 



W(L,t) 




(35) 



From equilibrium phase transitions, say in a magnetic 
system, we know that the fluctuations in the overall mag- 
netization is proportional to the susceptibility, which is a 
response function. It is known that this quantity diverges 
at the critical point for second order phase transitions, 
while it stays finite for first order phase transitions. For 
a d-dimensional system of size L, the divergence appears 
as x{T c ) ~ M 2 /L d ~ L 2 ~ n , with r\ being a critical expo- 
nent. 

We can study the overall fluctuations in the order pa- 
rameter field for the coating transition 56(x, t) = 0(x, t) — 
9 in the contact line problem. One can define 

^ = LJ d ^ dV (50(x, t)S9(x',t)), (36) 

and use the relation between 9(x,t) and h(x,t) given in 
Eq. ©, to find 
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j k 2 e? k ~ (\h(k,t)n (37) 



(27T) 



to the leading order. We can now use the scaling form of 
Eq. (|33|) in the long time limit, and obtain 
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which yields 
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We thus find interestingly that the value of the rough- 
ness exponent at the onset of the coating transition can 
determine the order of the dynamical phase transition: 
the phase transition is first order for £ < 2 -^ L : and it is 
second order for £ > 2 -^. In the case of a second order 
phase transition, we can define a dynamical exponent 

r} = 4 - d - 2C, (40) 

based on the above analogy with the equilibrium critical 
phenomena, and an order parameter exponent 



(41) 



from a naive application of the exponent identities, where 
v is the exponent that characterizes the divergence of the 
correlation length (see below). 

Similarly, we can calculate the magnitude of local order 
parameter fluctuations. We find 



to the leading order, which yields 

dk 



(59(x,t) 2 
and, consequently, 

<<W(x,i) 2 

where 



n/L 



k 2 i 



^9(kH). 




(42) 



(43) 



(44) 



r 1 



69 st (L) ~ { y/]n(L/a) ; ( = 1 



(45) 



C>1 



The above result shows that the extent of the order pa- 
rameter fluctuations actually depends on the value of the 
roughness exponent £ in the stationary limit, i.e., it will 
be finite for £ < 1, and unbounded for £ > 1. It is im- 
portant to note that £ > 1 signals a breakdown of our 
perturbative scheme in dealing with the nonlinearities in 
the system, as one can show that the neglected nonlinear 
terms in Eq. (|2^) can actually be organized as a power 
series in the parameter L'* -1 in the long wavelength limit. 

A systematic study of the stochastic dynamics de- 
scribed by Eq. (H|) will yield the values of z and C, and 
thus provide us with a characterization of the statistical 
properties of the moving contact line. 
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V. CONTACT-ANGLE-VELOCITY RELATION: 
MEAN-FIELD THEORY FOR THE COATING 
TRANSITION 



III , can be first studied 



The dynamics described in Sec. 
in the mean-held approximation, where the contact line 
is assumed to be a straight line. In this case, the relation 
between the dynamic contact-angle of the wetting front, 
and its velocity is given by Eq. (|25|). It is instructive to 
examine the limiting behavior of this equation in the two 
cases of local dissipation and hydrodynamic dissipation 
scenarios, separately. 



A. Local Approach 



In this regime, Eq. ( |25| ) yields [see also Eq.(Bl)] 



for v < £ j-, while 8 



co/ 

identically for v > 



(46) 



This 



function is plotted in Fig. ||. 

As can be readily seen from Fig. 0, increasing v would 
lead to decreasing values of 8 until at a critical velocity 
v c i = -g 1 it finally vanishes continuously. A vanishing 
contact-angle presumably corresponds to formation of a 
Landau-Levich film. The value of the dynamic contact- 
angle 8 serves as the order parameter for this dynamical 
phase transition, while v is the tuning parameter. The 
continuous vanishing of the order parameter makes the 
phase transition classified as second order. As in the gen- 
eral theory of critical phenomena, a mean-held exponent 
(3 = 1/2 is characterizing the vanishing of the order pa- 
rameter in terms of the tuning parameter. 



(9/9) 




FIG. 5. The reduced order parameter (6/6 e )i as a function 
of the dimensionless velocity v/coi for local mechanism [Eq. 
(g^)]. The dynamical phase transition at v c i/cot = 1/2 is 
predicted to be of second order in this picture. 



B. Hydrodynamic Approach 



In this case, Eq.fl25|) leads to [see also Eq. (B7)] 1 

1/3 




(47) 




FIG. 6. The reduced order parameter (0/6 e )h as a func- 
tion of the dimensionless velocity v/coh for the hydrodynamic 
mechanism [Eq. (p7|)]. The dynamical phase transition at 
Vch/coh = l/(3\/3) — 0.192 is predicted to be of first order in 
this picture. 

The above formula, which holds only for v < has 
two branches and only the one that recovers 8 = 8 e for 
zero velocity is acceptable as plotted in Fig. ^. While at 
v = ■jj s % we find = , we expect to have 8 = for 

3v 3 v 3 

higher velocities v > |2t. Therefore, the order parame- 
ter 8 experiences a finite jump at the transition velocity 
v ch = which is the hallmark of a first order phase 
transition. 



VI. LINEAR THEORY 

One can go beyond the simple mean-field treatment, 
and study the effect of linear fluctuations on top of the 



x Note that the expression in Eq. (^t|) is real, and the i is 
retained only to keep the appearance of the formula simpler. 
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mean-field theory. There are different aspects to the lin- 
ear dynamics that one can elaborate on, as considered in 
this section. 



A. Relaxation 

We can study the relaxation dynamics of a moving con- 
tact line only by using the linear term in Eq. (|23|). We 
find that each deformation mode of wavevector k relaxes 
with a characteristic time scale given as 



C. Breakdown 



- 1 (k) = c\k\, 



(48) 



where the relaxation vel ocit y c is give n as in Eq. ( |2q ) 
above. [See also Eqs. (B2) and (B8) for the limiting 
forms.] 

In fact, one can show from Eqs. ( |25| ) and (|26|), that 
the relaxation velocity c is related to the slope of the 
contact-angle-velocity curve as 



(49) 



This proves that the onset of instability in the contact- 
angle, signalled by a diverging slope of d9/dv, exactly 
corresponds to where the relaxation velocity c vanishes. 
In other words, exactly at the onset of the coating transi- 
tion (as found in the mean- field scheme of Sec. [v]), where 
the order parameter has a singular behavior, the relax- 
ation becomes infinitely slow, and a distorted the contact 
line does not relax. 

This "coincidence" suggests strongly that the coating 
transition — the onset of leaving a Landau-Levich film — 
can be described as a dynamical phase transition in terms 
of the deformations of the contact line and its statistical 
roughness on a disordered substrate. 



B. Fluctuations 

To account for the statistical fluctuations of the mov- 
ing contact line on a disordered substrate, we can also 
include the noise term in Eq. (|23|), and calculate the 
width of the contact line. We find 



y/l, 



t < * 



W(L,t) 



ln[f], §<t<f, 

L 



(50) 



t > 



There are two important time scales in the above equa- 
tion: (i) the "microscopic" time a/c that corresponds to 
the crossover from local diffusive dynamics to collective 
motion along the contact line, and (ii) the "macroscopic" 
time L j c that corresponds to the crossover to the station- 
ary state. We thus find from Eq. ([30j) that the stochas- 
tic deformations of the contact line are characterized by 
Co = and zq — 1, within the linear theory. 



The nonlinear term in Eq. ( |23| ) will modify the above 
results if it becomes appreciable at long length scales, as 
compared to the linear term. To examine under what 
conditions this may take place, we can estimate the ratio 
of the two terms in Eq. (23t), which scales as 



c a (L/a 



,2Co-l 



aco^/ln (L/a) 



c(L/a)to 



Lc 



where c could be set by either cqi or cq^- The non- 
linear term is thus appreciable only when the smallest 
time scale in the linear theory a/c becomes comparable 
to L/co, which happens nea r the coating transition when 
c becomes small (see Sec. VIA ). This confirms that a 
consistent study of the coating transition should take a 
proper account of the nonlinear terms. 



VII. EFFECT OF THE NONLINEAR TERM 

Let us now attempt to systematically study the dy- 
namical phase transition in the moving contact line using 
an RG scheme. The dynamical equation, which can be 



generally written in d dimensions as 1 33 
d t h{k, t) = -ck h(k, t) + rj(k, t) 



y^A(q > k-q)Mq ) *)Mk-q,*), (51) 



with A(q, k — q) given as in Eq. (|24|), belongs to the 
general class of Kardar-Parisi-Zhang (KPZ) equations 
pi] p9| . We take a noise spectrum given by 

MM)) =o, 

(f/(k, t)r)(k', t')) = 2D(2%) d 6 d (k + k')5(t - t'), (52) 

and employ standard RG techniques following Ref. |3f| 
to calculate the RG equations describing the flow of the 
coupling constants. 



A. Perturbation Theory 

To perform this calculation, we first need to construct 
a perturbation theory that takes into account the non- 
linear term. This is done more easily if we perform a 
Fourier transformation in the time variable in Eq. (|5T 
which yields 



iuh{k, ui) 



—ck h(k, u>) + 7/(k, ui 
1 f dil d d q 



A(q,k- q) 



2?r (27r) d 
x /i(q, Q)h(k — q, u — O). 



(53) 
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This equation can then be re-written in the form 

h(k,u) = G (k,w)r?(k,o;) 

dQ d d q 



1 , N r dn d d q_ , 



2?r (2tt) 

x ft(q,fi)ft(k-q,o;-fi), (54) 
in which the bare response function is given as 

G (k,Lu) = — — . (55) 



ck — iuj 

We can then define the full response function G(k, ui) via 

ft,(k,w) =G(k,a;)r?(k,w), (56) 

and use the spectrum for the Fourier transform of the 
noise as 



(r?(k,a;)) = 0, 

(r?(k, w)»/(k / , w')) = 2L>(27r)' i+ V(k + k')<5(w + a/), (57) 

to find the response function perturbatively in the A- 
parameters. 

Up to second order in the perturbation theory, we find 



G(k, w) = G (k, ui) + 4 x -G (k, uj) 2 x 2D 
fdSl d d q 

x G (q, n)G (-q, -n)Go(k - q, w - 0), (58) 



which can be re-written as 
d9sm d ~ 2 9 



c^/q 2 + k 2 - 2kqcos9 - i{u - Q) (c 2 q 2 + fl 2 ) 

Xi(q 2 - qkcos9) + X 2 q\/q 2 + k 2 - 2kqcos9 + X 3 kq + A 3 fc\/g 2 + k 2 — 2kqcos9 - X 3 k 2 



X 1 kqcos9 + X 2 qk + X 3 q\/q 2 + k 2 - 2kqcos9 + X 3 k^ ' q 2 + k 2 - 2kqcos9 - X 3 (q 2 + k 2 - 2kqcos9) , (59) 



where A = tt / a is an upper cutoff for the wavevector set 
by an inverse microscopic length scale, and Sd is the sur- 
face area of a unit sphere in d dimensions. We can then 
perform the frequency integration, expand the integrand 
in powers of k/q (that is justified because we are inter- 
ested in the long wavelength behavior of the system), and 
integrate over the angular variable, to obtain 



G'\Klo) =Go 1 (k,u,)- 



D Sa 



2c 2 (2n) d 

x (Ai + A 2 )(A 2 + A 3 ) (J q d dq^j k. (60) 

Assuming a same form for the full response function as 
the bare one, namely 



G(k,w) 



1 



CRk — iuj ' 

one then obtain a renormalizcd clastic constant as 



(61) 



cr = c 1 1 - 



(Ai+A 2 )(A 2 + A 3 )J 
2c 3 



S d 
(2ir) d 



q d dq 



(62) 



A similar calculation can be performed to obtain the 
renormalized noise amplitude, that is defined via 

(h(k, uj)h(-k, -u)) = 2D R G(k, u)G(-k, -uj). (63) 

one obtains 



2D R = 2D + 2x x (2D) 2 J 



dCl d d q_ 



2ir (2ir) d 

x A(q, k - q)A(-q, q - k)G (q, O)G (-q, -SI) 
X G Q (k-q ! w-O)G (q-k,n-w) ! (64) 



which in the small k limit yields 



D R = Dll 



(Ai 



X 7 ) 2 D 



4c 3 



S d 
(2tt) c 



q d dq 



(65) 



Finally, one can show that similar to original KPZ 
problem |34| , none of the A-parameters are renormalized, 
so that we have 



A H (q,k-q) = A(q,k-q). 



(66) 



We can now use the results of the perturbation theory, 
to construct a perturbative RG scheme. 



B. Renormalization Group Calculation 

In order to recapitulate the perturbation theory into 
an RG calculation, we should only integrate out over 
a layer of wavevectors from A/6 to A an see how the 
coupling constants are affected by that. This step, 
which makes up the coarse graining procedure, will lead 
to the same results as in Eqs. (B2h, (pa), and (pfl), 



in which the wavevector integration reads J A y b q d dq 
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(l - This should then be followed by the 

scale transformations x — > bx, t — ► b z t, and h(x,t) — > 
b^h(x,t), which for the scale factor of the form b = e l 
and for infinitesimal values of / yields the following RG 
flow equations for the coupling constant: 



dc 
dl 

dA(q, k q 



= c[z-l-U], 



dD 
~~dT 

in which 



dl 
= D 



U 



2(-d 



-q) (C + ^ 


-2) 


/Ai + A 2 \ 


U~ 


VA2 + A3J 


~2 



(67) 



7r5 d (Ai+A 2 )(A 2 +A 3 )L> 



(2a) 



d+l c 3 



(68) 



To study the fixed point structure of the above set of 
flow equations, we set z = 1 + U and £ = 1 — U, and look 
at the flow equation for U: 



-=-(d + l)U + 



6 



Ai + A 2 
A 2 + Aa 



U 2 

T' (69) 



which has two stable fixed points at f7 = (linear theory) 
and U = 00 (strong coupling) , as well as an intermediate 
unstable fixed point at 



U = U* 



2(d + l) 



6+(A!+A 2 )/(A 2 + A 3 ) 



(70) 



For U < U* , the nonlinearity is irrelevant and the ex- 
ponents are given by the linear theory, i.e. £0 = and 
z = 1, while for U > U* the behavior of the system is 
governed by a strong coupling fixed point which cannot 
be studied perturbatively. The fixed point at U* corre- 
sponds to a roughening transition of the moving contact 
line. The exponents at the transition are 



z = i + ir, 
C = i-tr, 



(71) 



which turn out to be nonuniversal. The roughening tran- 
sition corresponds to the limit of stability of the moving 
contact line phase. The phase described by the strong 



coupling fixed point could presumably correspond to a 
Landau-Levich film or a pinned contact line. 

We can also study how the transition is approached by 
linearizing the flow equation near the fixed point. Setting 
U = U* + SU, wc find dSU/dl = (d + 1)SU that would 
imply divergence of the correlation length near the tran- 
sition as 



£~\5U\-», 



with 



d+1 



(72) 



(73) 



The correlation length corresponds to the typical size of 
rough segments in the contact line, which should diverge 
as the transition is approached (See Fig. Fl). 



VIII. FIXED POINT EQUATION: PHASE 
DIAGRAM AND EXPONENTS 

In principle, the position of a phase boundary that sep- 
arates the different phases could be obtained from a fixed 
point equation [such as Eq. (]70|)]. However, in most RG 
studies the relation between the phenomenological pa- 
rameters in the theory and the microscopic parameters 
are not known, and thus the fixed point equation cannot 
help us obtain the phase diagram of the system in terms 
of the real control parameters. 

In the present case, however, the fact that we have 
used some physical models to arrive at the the dynami- 
cal equations allows us to make such direct connections. 
If we use the relations obtained for the parameters as a 
function of the contact-angle in Sec. Ill (and Appendix 
|b|), and insert them in the fixed point equation Eq. (|70| ) 
with d = 1, which can be written as 



D 
~7a 



8 



7T (Ai + A 2 )(Ai +7A 2 +6A3) 



(74) 



we can map out the phase diagram of the system, and 
calculate the value of the exponents. 

In the local case, using Eqs. (Bl-B6)inEq. (^3) yields 



g 



(f) 1/2 (£ 



an equation for the phase boundary as 

1/2 



1 — 2- 



1 - (1 - 2fei + 2& 2 + 26g)f| 1 - (1 - 2&i + 146; 



1/2 



Ubl 



126 3 )fl 



1/2' 



(75) 



Similarly, for the hydrodynamic case, using Eqs. (B7-B12) in Eq. ( |74| ) yields an equation for the phase boundary as 

3/2 



e h 



(i 



(?) 



3 X l/2 



0?. 



1/2 



(61 - 6a - 1) + (1 - 36i + 36a + 6fog)|^ _ 7b 2 - 66 3 - 1) + (1 - 36i + 216 2 + 426jj + 1863)^ 



1/2 r 



1/2 ' 



(76) 
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These phase boundaries are plotted in Figs. ^ and || (solid line) for a choice of parameters bo = 1, 61 = 1, 62 = — 1, 
and 63 = — 1. 

The phase boundary that corresponds to a roughening transition of the moving contact line has two different 
branches. The first branch that happens at relatively high velocities, presumably corresponds to the onset of leaving 
a Landau-Levich film. In the local case, the phase boundary for this transition start at zero contact-angle for zero 
disorder, and has a limiting form as reported in Eq. (Ill) above, with 



Oicl 



(27T) 



1/6 



(77) 



whereas in the hydrodynamic case the boundary starts at a finite value of the contact-angle, with the limiting form 
as reported in Eq. (0), in which 



1 



OLch 



2V3&o 



-(3&g-l)(216§-l) 



1/3 



(78) 



We can calculate the value of the exponents along these phase boundaries, and as it turns out they are non-universal. 
We find 



^ = 1 + 2(2^ ('£±»L 



: cl = 1-2(2^/3 f'M+ & 2 + &3 



b 



9 

9 

181 



2/3 



2/3 



for the local case, and 

3-1/(362) (361T) 1 / 3 (362-l)V3 



Cch 



Zch 



7-l/(36g) bo (2162-1)5/3 

11 - l/(3b§) (367^/3 (36^1)1/3 
7-l/(36g) + 6 (2162 - 1)5/3 



[362(6 1 + 6 3 -l)- 6 2 - 6 3 ] 
[36^(6i + 63 - 1) - 6 2 - 6 3 ; 



9 

9 



2/3 



2/3 



(79) 
(80) 

(81) 
(82) 



for the hydrodynamic case. 

The other branch of the phase boundaries appear very 
near the equilibrium contact-angle corresponding to low 
velocities, and presumably corresponds to the onset of 
pinning. In other words, the pinning transition (as op- 
posed to the depinning transition) can be described as a 
roughening of the contact line imposed by the static min- 
imum energy configuration on the disordered substrate, 
as it slows down to rest. The shape of the phase bound- 
aries for weak disorder can be found for the local case as 
in Eq. (||) above, with 



3?r /i 2 7 



1 /6n + fo 2 + 6 3 



6 V&0 + 62-61 



(83) 



and, for the hydrodynamic case as in Eq. (0) above, with 



37T 

a p h = ^3(^0 + 62-61) 



362 



6 2 + 63 



(84) 



The above prediction for the phase boundary of the pin- 
ning transition is in agreement with a previous predic- 
tion based on the hysteresis in receding and advancing 



contact-angles |iq] . We can calculate the corresponding- 
exponents along the depinning transition phase bound- 
ary, and find 



1 



(pi — 



Cph 



1 (bo 

2 U 



+ &2-f>l \ 



+b 2 +b 3 ) 



2 \bl+b 2 +b 3 



a , 1 ( tf,+b 3 -bi \ ' 
J+ 2 \bl+b 2 +b 3 J 



for the local case, and 



1 / 3fcg+fc 2 -bl 

2 1 36* +62+63 



1 (3bl 
2 



+6 2 -6i 
+b 2 +b 3 



3feg+6 2 -6i 
36j* +62+63 



1 (3bl 

2 {sbj, 



+6 2 -&i 



+62+63 



for the hydrodynamic case. 



(85) 



(86) 



(87) 



(88) 
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The two branches of the phase boundaries meet at a 
junction point, called T in Figs. || and |[ where they both 
develop a vanishing slope as obtained from Eqs. ( |r5| ) and 

Although the fact that we do not know the values of 
the numerical constants b n makes us unable to predict 
the values of the exponents, we can nevertheless put some 
bounds on them using simple physical requirements. For 
example, we expect on physical grounds that the critical 
contact-angle for the coating transition increases with 
disorder. Then Eqs. (||) and (|78|), together with a re- 
quirement that we have a roughness exponent that is 
less than one, imply that fe 2 , > 1/3, which results in a 
criterion 



1 /• 3 

g < (,ch < 



(89) 



for the roughness exponent in the hydrodynamic case. In 
the local case, to have a roughness exponent that is less 
than one we should have 6q + &a + &3 < 0, as can be seen 
fromEq. @. 

From the above results and bounds on the values of 
the roughness exponents in the two different approaches, 
and the analysis of Sec. IV, we can conclude that at 
least in the weak disorder limit £ c ; > | and Cch < \ , and 
thus the coating transition remains to be first order in 
the hydrodynamic picture and second order in the local 
picture even beyond mean-field theory, i.e. within the 
RG scheme. 

On the pinning transition boundary a similar require- 
ment that the advancing (receding) contact-angle in- 
creases (decreases) with disorder (together with a re- 
quirement that we have a roughness exponent that is less 

than one) leads to ( l^^ 2 , f 1 ) > 0, as can be seen from 

\ Og+02+03 J 



Eqs. (H) and (183) , which implies 



1 



(90) 



and ^ 3^^ + b3 ) > 0, as can be seen from Eqs. (||) and 



^3bl+b 2 +b 3 J 

(|S4|), which implies 



Cph > -■ 



(91) 



The above results are in agreement with the experimen- 
tally observed roughness exponent of about 0.5 for both 
liquid Helium on a Cesium substrate |25| and water on 
glass experiments pq] . 



IX. DISCUSSION 

The above study of the dynamics of a moving contact 
line on a disordered substrate reveals that both the in- 
stability corresponding to the onset of leaving a Landau- 
Levich film at high velocities, and the pinning of the 



contact line to the substrate at low velocities, can be de- 
scribed in a unified framework as roughening transitions 
in the contact line. For advancing contact lines, the phase 
boundary corresponding to pinning extends continuously 
to the strong disorder regime. For the receding case, how- 
ever, the phase boundaries corresponding to pinning and 
to Landau-Levich transition asymptote to a maximum 
with zero slope, where we identify as a junction point, 
denoted by T in Figs. ^ and |[ 

For stronger disorder, the pinned contact line presum- 
ably goes directly into the Landau-Levich phase: the 
edge of the liquid drop is pinned so strongly that it re- 
mains still when the liquid wedge starts to move upon 
decreasing the contact-angle, and thus a film is left be- 
hind. In this sense, at the dashed line in the phase 
diagram (Figs. |2| and ||) nothing really happens to the 
contact line itself, while the body of the liquid drop is 
"depinned." This behavior has in fact been observed in 
experiments with liquid Helium on a Cesium substrate, 
where the receding contact-angle has always been found 
to vanish identically immediately after depinning p8| , p5[ . 
It is also interesting to note that a non-universality of the 
roughness exponent at the depinning threshold, manifest- 
ing itself in the form of a dependence on the value of the 
contact- angle, has also been reported in Ref. [^5| . 

One can get a plausible picture of the roughening tran- 
sition in terms of fluctuating domains of different sizes, 
in analogy with equilibrium phase separations in, say, bi- 
nary mixtures. As it is sketched in Fig. a portion of 
the moving contact line can be instantaneously pinned to 
a defect on the substrate, thereby nucleating a domain 
where the liquid is pinned to the substrate. These do- 
mains are rough, because the contact line has to conform 
to the minimal energy configuration imposed by the sub- 
strate disorder, and they can exist in length scales up to 
the correlation length £. As the transition is approached, 
the correlation length increases until at some point one 
of the domains enlarge to a macroscopic scale and span 
the whole system, corresponding to the divergence of the 
correlation length at the critical point. Interestingly, the 
same picture can be attributed to both the coating tran- 
sition and the pinning transition. The difference comes 
from the fact that in the pinning transition the substrate 
can stop both the contact line and the liquid from mo- 
tion, while in the coating transition where the liquid is 
already moving at relatively high velocities it can only 
stop the contact line and not the liquid, as mentioned 
above. 

One can use a scaling argument to account for the 
power law in Eqs. (§) and (§). If we take the above KPZ 
equation and make the scale transformations t — > ct, 
and h — > ^D/ch, the coefficient of the linear term as 
well as the strength of the noise term will be set to one, 
and the only remaining coupling constant (the coefficient 
of the nonlinear term) will have the form AD 1 / 2 /c 3 / 2 . 
We expect the roughening transition to take place when 
this coupling constant is of order unity, which yields 
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Eqs. (|T]) and (Q) when the values D ~ g 2 , c; ~ |§, 




and A ~ const, are used. 



In treating the disorder contribution, we have made the 
assumption that the dependence in the noise term r\ on 
the shape of the contact line can be neglected for any non- 
vanishing drift velocity. While this approximation can be 
supported by a naive power counting, it may be argued 
that more a sophisticated functional RG approach is nec- 
essary to deal with this dependence [ pl"|]23| , since upon 
approaching the pinning transition the contact line veloc- 
ity becomes progressively small. In this sense, the exact 
values of the exponents derived for the pinning transi- 
tion boundary, which are non-universal anyway, and the 
predictions about the nature of the junction point may 
not be trusted. However, we do believe that the main 
features of the above results — the non-universality of the 
exponents, the existence of the junction point, and the 
form of the pinning phase boundary, will persist. 

We finally mention that this work can hopefully mo- 
tivate two types of experiments. One can try to probe 
the relaxation of a moving contact line, similar to the 



Ondarcuhu-Veyssie experiment on a static contact line 
[fiJI, and measure the velocity dependence of the dis- 
persion relation. This dependence could be used to de- 
termine the dominant dissipation mechanism p^j . It is 
also interesting to systematically study the coating tran- 
sition (onset of leaving a Landau-Levich film) for reced- 
ing contact lines on a disordered substrate using video 
microscopy techniques. In particular, it would be inter- 
esting to look for a roughening of the contact line before 
the Landau-Levich film is formed, and measure the cor- 
responding roughness exponents. 
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APPENDIX A: GEOMETRY OF THE SURFACE NEAR THE CONTACT LINE 

Here we briefly review some of the aspects of the differential geometry of surfaces that are useful in this context. 
The free surface of the liquid can be described by the embedding 

R(x, y) = (x,y,z(x,y)), (Al) 

using which we can find two independent unit tangent vectors to the surface at each point 

' " 1 =(M, (A2) 



\o x n\ ^1 + (d x z)- 

and 



^^TttW™" (A3) 



We can also define the unit tangent vector for the contact line (see Fig. ^) as 

1 



y/1 + (W 



(l,cU,0). (A4) 



To find the unit vector T, that shows the local direction at which the surface tension force is acting (see Fig. ||), we 
should note that it is tangent to the surface, so it can be written as 

T = u x i x + Uyiy. (A5) 

Imposing the requirements that it is perpendicular to the contact line T • t = and that it is a unit vector T 2 = 1, 
yields the two parameters, and we obtain 

rf, _ (-d x h,l,-d x hd x z + d y z) 



{i + {d y zf + {d x hy [i + {d x zY\ - 2(d x h)(d x z)(d y z)} 1/2 

We can also define the unit vector normal to the contact line (see Fig. 0) as 
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1 



n 



= (-cU,l,0)), 



using which we can define the local contact-angle as 

6{x) = cos' 1 (± ■ nj , 

where 



T n 



y/1 + (d x hy 



{i + {d yZ y + {d x hy [i + {d xZ y\ - 2{d x h){d x z){d yZ )} 



1/2- 



(A7) 



(A8) 



(A9) 



Note that we have v/l + (d x h) 2 T ■ y = cos 6(x) 



APPENDIX B: COUPLING CONSTANTS IN THE 
TWO DIFFERENT APPROACHES 



2. Hydrodynamic Approach 



The corresponding asymptotic forms of the equations 
In this Appendix, we have summarized the limiting giyen in gec g in thig cage can be obtained by taking 



form of the coupling constants of the dynamical equa- the Hmit Coh/m ^ Q Qne obtaing 
tion, corresponding to each of the dissipation mechanisms 
separately. 



coh 



1. Local Approach 

The asymptotic form of the equations given in Sec. 



[II can be obtained in this case by taking the limit 
coh/coi — > co. One obtains 



9 1 
0* 



for the contact-angle-velocity relation, 

Q = boCol—r, 

for the relaxation speed, and 



Aiz- T 



(Bl) 



(B2) 



(B3) 



and 



and 



A 2 ; = (b 2 + b )c i-^, 



^31 = b^CQi—r, 



(B5) 



for the nonlinear coupling constants. We can also find 
the strength of the noise term as 



«=— (4)' 



(B6) 



for the contact-angle-velocity relation, 



c h = bo — j- (3^-1 



for the relaxation speed, and 



(6i-l) + (1-360- 



and 



and 



2 fl 



-6a + (36 2 + 66g) 



q2l 



^ 2 



A 3 ft = b 3^-j- ( - 1 S • 



(B7) 



(B8) 



(B9) 



(BIO) 



(Bll) 



(B4) f° r the nonlinear coupling constants. We can also find 
the strength of the noise term as 



D h = 



a 2 c h ( 6 



1 - 



9 

101 



(B12) 
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